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Abstract 



The optical model potentials for nucleon-nucleus elastic scattering at 65 MeV 
are calculated for ^^C, ^^O, ^^Si, ^°Ca, ^^Fe, ^OZr and ^o^Pb in first order mul- 
tiple scattering theory, following the prescription of the spectator expansion, 
where the only inputs are the free NN potentials, the nuclear densities and 
the nuclear mean field as derived from microscopic nuclear structure calcu- 
lations. These potentials are used to predict differential cross sections, ana- 
lyzing powers and spin rotation functions for neutron and proton scattering 
at 65 MeV projectile energy and compared with available experimental data. 
The theoretical curves are in surprisingly good agreement with the data. The 
modification of the propagator due to the coupling of the struck nucleon to 
the residual nucleus is seen to be significant at this energy and invariably 
improves the congruence of theoretical prediction and measurement. 
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I. INTRODUCTION 



Recently measurements of 65 MeV neutron-nucleus elastic differential cross sec- 
tions have been published Together with the corresponding proton data, one now has 
an opportunity to study the effects of isospin degrees of freedom and to analyze multiple 
scattering theory with and without the coulomb field at low energies. In the past it had 
been assumed that first order multiple scattering theory would prove unable to provide an 
accurate representation of experiments at such low energies and thus models in the form of 
phenomenological optical potentials or effective nucleon-nucleon interactions were proposed. 
With recent advances in our ability to calculate the first term in a multiple scattering expan- 
sion together with the the influence of the mean field potential binding the struck nucleon 
to the target nucleus, it is now possible to test whether one is able to represent low-energy 
scattering in the first order or whether one must consider higher order scattering terms, all 
within a parameter-free description. 



II. FORMALISM 

We have been pursuing a program of calculation of elastic nucleon-nucleus scat- 
tering at energies sufficiently high such that first order multiple scattering theory in the 
forward cone provides a good description of the data [0,0. In this program the only inputs 
are the nucleon-nucleon (NN) interaction, represented by NN t-matrices, the target wave 
functions, and the static target nuclear mean field. These quantities are incorporated into 
multiple scattering theory in the hierarchical spectator expansion for the optical potential, 
in which the transition operator, T, is defined to be 

T = U + UGoPT. (2.1) 

The operator P is the projector onto the target initial state. The optical potential operator 
U is given in the spectator expansion as: 
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^ = E^^+ E^^. + --- • (2-2) 

The first order theory corresponds to a truncation of this series to a single term, 

f/^E^- (2-3) 

i 

where 

Ti = voi + VoiGo (1 - P) Ti (2.4) 

Here Voi stands for the potential between the projectile nucleon (0) and the zth target 
nucleon. The propagator Gq is given by 

Gq^ = E — ho — Ha 

= E -ho-hi-Y^^ij - 

= E-ho-H,-H\ (2.5) 

where hi is the kinetic energy operator for nucleon i and = Hi + is the target 
Hamiltonian. The operator, Tj, can be reexpressed in a solvable one-body integral equation 



as 



Ti — ioi ~ toiGoPTi , (2-6) 



where toi is given to be: 



^oi = toi + toi [Go — go\ toi 

= tlr^ + tiT^%%9otoi . (2.7) 

Here t^l^'^ is the free NN t-matrix and go is the free NN propagator for the active pair 
consisting of projectile and target nucleon g^^ = E' — h^ — hi . The scattering operator 
% expresses the scattering of the target nucleon (z) from the residual nucleus, which is 
represented by J^jjti'Vij. The explicit treatment of Eq. { \i.7\ ) is described in detail in Ref. |0] 
and is directly derivable within the spectator expansion of multiple scattering theory. The 



right term of Eq. ( p.7|) results from the difference between the free propagator go with Go, 
which corresponds to the propagation of the target nucleon through the nuclear medium 
and can be thought of as a propagator modification. 

The first order optical potential is then constructed with the operator, Tj, from Eq. ( ^.61 ): 

Uopt = {k'o^A\T.^A^A) . (2.8) 

i 

In the present calculations, which are performed in momentum space, Uopt enters in the 
'optimum factorized' or 'off-shell rp' form as 

Uopt^r{q,}C-E)p{q) , (2.9) 

where q = k'Q — ko and /C = | {k'^ + ko); k'^ and /cq are the final and initial momenta of the 
projectile. This corresponds to a steepest descent evaluation of the 'full-folding' integral, in 
which the r is convoluted with the nonlocal density as indicated schematically in Eq. ( p.SI) . 
For harmonic oscillator model densities it has been shown that the optimum factorized form 
represents the nonlocal character of Uopt qualitatively in the intermediate energy regime . 
Complete 'full-folding' calculations with more realistic nuclear densities are in progress. It is 
to be understood that all spin summations are performed in obtaining Uopt (under the usual 
assumption of a spin-saturated target), thus reducing the required NN t-matrix elements 
to the spin-independent component (corresponding to the Wolfenstein amplitude A) and 
the spin-orbit component (corresponding to the Wolfenstein amplitude C). All scattering 
calculations presented here contain an additional factor in the optical potential to account 
for the transformation of the NN t-matrix from the two-nucleon cm. frame to the nucleon- 
nucleus cm. frame 0. 

III. THEORETICAL PREDICTIONS 

For a calculation of the first order optical potential in the optimum factorized 
form, the quantities T{q,IC;E) and p(g) are required as input. All calculations presented 
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in this paper are based on the full Bonn Potential as the NN interaction from which 
r(g, /C; E) is obtained. As can be seen from Eq. ( |2.6| ), the quantity (rp) can be calculated as 
the solution of a one-body integral equation in which {tp) serves as the driving term. The 
proton densities are taken from charge densities measured in electron scattering experiments 
Although the neutron distribution in nuclei are not completely determined by measure- 



ment, nuclear structure calculations indicate significant differences between the proton and 
neutron distributions. These differences are nonnegligible in our reaction calculations and 
especially affect the location of the diffraction minima and become even more pronounced in 
the case of neutron scattering. We believe neutron densities taken from the Hartree-Fock- 
Bogolyubov calculations of Ref. [§] give the best presently available representation of the 
neutron distributions. These are the densities used in the present work. 

The calculation of toi according to Eq. ( p.7|) requires the free NN t-matrix as well as 
%, the 't-matrix' representing the scattering of the struck target nucleon from the residual 
nucleus. A one-body mean field potential for a nucleon within the target nucleus is used 
as the driving term to obtain % as solution of a Lippmann-Schwinger type equation. The 
inclusion of this correction, 7^, corresponds to a modification of the free NN propagator, qq, 
due to the nuclear medium in order to recover the propagator Gq, which correctly represents 
the propagation of the target nucleon through the nucleus. Our calculations use two different 
models for the mean fields, one is the nonlocal, nonrelativistic mean field potential taken 
from a Hartree-Fock-Bogolyubov calculation IQ . Curves based on this choice are represented 
by dash-dotted lines. The second choice involves a nonrelativistic reduction of the mean field 
potentials resulting from a Dirac-Hartree calculation based upon the a-u model [|10| . Curves 



based on this choice are represented as dashed lines. Calculations using only the free NN 
t-matrix correspond to a truncation of Eq. ( p.7| ) after the first term, so that 

io^ ^ ir ■ (3.1) 

Results using this truncation are included in the figures to illustrate the importance of the 
correction due to % in the energy regime under discussion. The corresponding curves are 
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shown as solid lines. 



IV. COMPARISON WITH DATA 

The elastic neutron and proton elastic scattering observables [differential cross 
section analyzing power Ay and spin rotation function Q] are calculated for ^^C, ^^O, 
28gj^ ^°Ca, ^^Fe, ^^Zt and ^''^Pb. In all of these cases, except for ^^O and ^°Zr, neutron 
differential cross section data ]l| exist for the natural elements along with the more extensive 
proton data. The carbon calculations in Fig. 1 should not be taken as seriously as the others, 
since it is known that carbon is a highly deformed nucleus. For a realistic calculation of 
carbon certain additional collective degrees of freedom need to be considered. 

The overall impression is that the first order, parameter free multiple scattering pre- 
dictions provide a very good representation of the data for angles below ~ 60°, which is a 
very pleasant surprise. One possible partial explanation for this good agreement is that the 
momentum transfer for 65 MeV for 9 = 180° is roughly equivalent to that corresponding 
to 9 = 60° at 260 MeV. That is to say, the range of momentum transfer under observation 
is rather small relative to that considered at higher energies. Nonetheless, this good a de- 
scription of the data is unexpected. At this low an energy, the higher order terms in the 
spectator expansion were expected to become important. It now appears that these terms 
will not be strongly structured and so will, no doubt, raise the forward cross section and fill 
in the sharp diffraction minima. A change in the spin structure due to higher order terms 
may only be apparent at large scattering angles. 

In a recent publication the authors concentrated on the total cross section atot for neutron 
scattering from ^^O and ^''Ca as a function of scattering energy It was found that our 
calculations reproduced the neutron total cross section data above ~ 100 MeV, whereas at 
65 MeV the first order multiple scattering theory predictions for the total cross sections for 
^°Ca fall short of the measured values by about ~ 10%. Thus, for example, if one looked 
at the proton elastic differential cross section in the forward cone, one might be surprised 
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to see how well the data is represented by the theoretical predictions. On the other hand, 
it is observed that the theoretical curve for the neutron elastic differential cross section 
consistently falls below the data in the forward cone. The reason is that in the proton case 
the Coulomb effect dominates at small scattering angles and thus masks any underprediction. 
It may also be argued that the neutron projectile should be able to penetrate further into the 
nucleus than a proton projectile, due to the lack of the coulomb barrier. This may perhaps 
mean that for heavier nuclei, neutron scattering might be more sensitive to higher order 
multiple scattering effects as well as to the medium corrections used here. This is offered to 
illustrate the enhanced value and impact of having combined neutron and proton data sets. 

The curves presented show very clearly that the inclusion of the propagator modifica- 
tion arising from the coupling of the struck nucleon to the residual nucleus improves the 
description of the data, in some cases dramatically. For neutron scattering only differential 
cross section data for angles up to 60 degrees are at present available. While this data is 
limited, a definite improvement from the propagator modification is seen. Relative to the 
unmodified calculation the diffraction minima systematically move to higher angles in bet- 
ter accord with the observation and the diffraction minima become less sharp also in accord 
with observation. For proton scattering, extensive data are available for Ay and Q. 
While there is always marked improvement in the differential cross section description due 
to the propagator medium modification, the effect on the spin observables Ay and Q is even 
more striking. For Ay at the first diffraction minima the propagator modification moves the 
theoretical curve closer towards the data so that for up to about 50 — 60 degrees there is a 
pretty good representation of the data. The effect on the Q function is even more dramatic 
where the medium correction shifts the curve down significantly so that the theoretical curve 
sits almost directly on top of the data for angles up to about 60 degrees. It appears that 
for Q, the propagator modification is required to obtain an excellent agreement with the 
data, and it is further remarkable that this is the case for all nuclei under consideration. 
The theoretical predictions for the heavier elements ^^Zi and ^'^^Pb are not as good as for 
the lighter nuclei, but the dramatic effect observed for the Q function remains. The effect 
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of the propagator modification on tlie spin rotation is slightly larger for neutron scattering 
than for protons and in Figs. 4-7 this difference is even larger for both Q and Ay. The good 
description of the spin observables independent of the mass of the nucleus may be taken 
as evidence that at these low energies the scattering process is generally completely surface 
dominated, and by construction our theoretical approach captures the effect of the coupling 
of the struck target nucleon to the mean field of the residual nucleus. 

All the nuclei under consideration are even-even nuclei. However, in a shell model sense 
they are mostly not 'spin-saturated', and so one might expect that all Wolfenstein parameters 
in the NN t-matrix to contribute to the optical potential. The remarkable description of 
the data which is obtained here based on only A and C suggests that such effects are small. 
This may give information concerning the structure of the target, but in the present stage 
it would be premature to test sensitivities to the correlation structure of the target. 



V. SUMMARY AND CONCLUSION 



The success of our calculations at higher energies [0,|ri|] have emboldened us to venture 
toward energies lower than previously considered. Neutron and proton elastic scattering 
observables are calculated in first order multiple scattering theory in a parameter free fashion 
for a number of even-even spin zero targets at 65 MeV. Even at this low energy a good 
description of the data is obtained. It is observed that the propagator modification due to 
the coupling of the struck nucleon to the residual nucleus invariably brings the theoretical 
results into closer agreement with observation. 

The availability of both neutron and proton scattering data at the same energy is es- 
pecially valuable. At this energy we find that our theoretical predictions consistently un- 
derpredict the neutron differential cross section in the forward direction; whereas, in the 
corresponding proton differential cross section. Coulomb effects mask this underprediction. 
This illustrates the enhanced value and impact of having available combined neutron and 
proton data sets. 
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FIGURES 

FIG. 1. The differential cross section, analyzing power and spin rotation functions are shown 
for elastic nucleon scattering from ^^C. The left three panels show elastic neutron scattering, while 
the right three panels are for proton scattering. The neutron data are from Ref. while the 
proton data are from Ref. In all cases, the solid line corresponds to a calculation of the 

first order optical potential based on the free NN t-matrix from the full Bonn model [Q] as input. 
The dashed and dash-dotted lines include the propagator modification due to the nuclear mean 
field. For the dash-dotted curve a Hartree-Fock-Bogolyubov ^ is used, for the dashed curve a 
Dirac-Hartree mean field |jlO|. 

FIG. 2. The same as Fig. 1, except for ^^O and the elastic proton scattering data are from 



Ref. 1141. 



FIG. 3. The same as Fig. 1 except for "Si and the elastic proton scattering data are from 



Ref. [|14|. 



FIG. 4. The same as Fig. 1 except for ^'^Ca and the elastic proton scattering data are from 
Ref. PI. 



FIG. 5. The same as Fig. 1 except for ^^Fe and the elastic proton scattering data are from 



Ref. [|14| 



FIG. 6. The same as Fig. 1 except for ^"Zv and the elastic proton scattering data are from 



Ref. 114 1 . 



FIG. 7. The same as Fig. 1 except for ^'^'^Pb and the elastic proton scattering data are from 
Ref. p|. 
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Figure 4 



Ca, neutron scattering Ca, proton scattering 
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Figure 7 



Pb, neutron scattering Pb, proton scattering 
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